Full two-dimensional numerical study of the quantum-Hall Skyrme crystal 
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Spin textures arise in the effective action approach to the quantum Hall effect. Up to now these 
textures, afso called Skyrmions, have been mainly studied using approximations. Studies of its fi- 
nite density form, the crystal, have been limited to a given symmetry and tend to ignore Skyrmion 
deformation. Using a simulated annealing technique, adapted to problems with a long-range inter- 
action, we are able to present the first full two-dimensional study of quantum Hall Skyrmions and 
the corresponding crystals. Our results show that only the Skyrmions with topological charge one 
and two are bound and that the Skyrme crystal is probably a triangular lattice made of charge-two 
Skyrmions. 
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In a strong magnetic field a two-dimensional electron 
gas exhibits all the complexity of the quantum Hall ef- 
fects 0. The motion of electrons is described by the 
Landau levels. At certain filling factors of the low- 
est Landau level (notably 1 and 1/3), the ground state 
becomes ferromagnetic ||. The lowest excitations are 
no longer single-particle states but spin textures j|, a 
type of solitons called Skyrmions. The Skyrmion is a 
twisted two-dimensional configuration of spins with spin 
down at the center and spin up far away from the cen- 
ter ||. A Skyrmion has a topological charge which 
corresponds to the number of times the spin rotates 
by 2ir. Its stability arises from the balance between 
the electron-electron Coulomb energy (of the order of 
e 2 /(47rd) w 13.7 meV and the Zeeman energy (of the 
order of gp B B w 0.3 meV) g, §. 

Since Skyrmions are the basic low-energy excitations 
in the quantum Hall effect at v = 1, they should play a 
crucial role in the groundstate near v = 1. It is expected 
that this state contains a finite density of Skyrmions || . 
Do these Skyrmions form a crystal? If so, what is the 
symmetry of the crystal? How important are quantum 
fluctuations and disorder caused by impurities? No fi- 
nal consensus as to the relative importance of all these 
effects has yet been reached. All experimental evidence 
is consistent with the existence of Skyrmions fit]. A re- 
cent paper by Khandelwal et al. [Bj gives evidence for 
Skyrmion localization, but does not seem to fully agree 
with most theoretical predictions || . 

One of the reasons may be that up to now the study 
of the Skyrme crystal and multi-Skyrmions, Skyrmions 
of arbitrary topological charge, has been performed us- 
ing simplifying approximations. Most calculations as- 
sume the existence of a symmetry and/or start from a 
simple potential picture that ignores deformation of the 
Skyrmions. For multi-Skyrmions, one uses an Ansatz. 
It is difficult to numerically minimize the Skyrmion en- 
ergy functional, because it contains a non-local Coulomb 
term. On the other hand, full two-dimensional nu- 
merical studies have been done for a local version of 



the two-dimensional Skyrme model, usually called the 
baby Skyrme model to distinguish it from the three- 
dimensional one. A thermodynamic approach was used 
to find its multi-Skyrmions (using simulated annealing) 
|]l0f and to study its thermodynamical properties [pd| . 
We have now successfully extended this promising tech- 
nique to non-local Coulomb terms. 

In this Letter, we are therefore able to numerically dis- 
cuss the structure of the quantum Hall multi-Skyrmions 
and Skyrme crystals without any restricting approxima- 
tion. We will also use our experience on the thermody- 
namics of the local baby Skyrme model to extrapolate to 
the quantum Hall model and discuss future research. We 
will give more details on numerical techniques and more 
detailed results in a longer paper Q. 

As mentioned, the quantum Hall ground state is fer- 
romagnetic for filling factors v = 1 and v = 1/3. 
The lowest Landau level projected quantum field theory 
is quite complicated and has interesting quasi-particle 
states. One can either try to use quantum-mechanical 
Hartree-Fock approximations 1 1 3 , or use a semiclassical 
Skyrme model. We use the latter approach, where the 
Skyrme model arises through the effective action of the 
long-wavelength limit of the quantum field theory. Its 
Lagrange density has a time dependent term containing 
a gauge interaction, a gradient expansion term, a Zeeman 
term, and a Coulomb interaction: 
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Here, Vym is the filling factor of the ferromagnetic compo- 
nent of the QH state. The ferromagnetic charge density 
is P — f^rl: where the magnetic length I = y/h/eB. The 
topological charge density, q(r) = n ■ (d x n x d y n)/ATT, 
describes the deviation of the electron density from the 
ferromagnetic background, p(r) — p = vpyiqif)- The spin 
stiffness p s = ,— vfma — r- The relative dielectric con- 
stant e r (e = e r eo) and the g factor are properties of the 
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sample used. 

In this Letter, we ignore the time-dependent term, be- 
cause we are looking for static solutions. Briefly, its 
role is to make the field n satisfy the Poisson bracket 
{rii(r), rij(r')} = Cijkrik(r)5(r — r'), and the Hamilto- 
nian is equal to the potential energy. 

The Zeeman term is proportional to the unit of length 
squared and the Coulomb term is proportional to the in- 
verse length d. The gradient term is scale invariant. The 
size of the static soliton solutions of the Skyrme model 
is thus determined by the balance between the Coulomb 
and the Zeeman term. Using a virial argument, one can 
show that at the minimum energy solution the Coulomb 
energy must be twice the Zeeman energy. We find it nu- 
merically preferable to use a unit of length in which the 
Zeeman and Coulomb terms balance, rather than use the 
magnetic length as unit. Equating the coefficients, we get 
gpfisBd 2 = Vp M e 2 / (Aned) , which leads to 
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They are our units of length, and the unit of energy be- 
comes e 2 /(47red). The energy density reads (x = r/d, 
B = Be 3 ) 
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This is thus effectively a one-parameter problem in terms 
°f Ps = 16 7g ^fm f ■ Note however that the quality of the 



gradient expansion depends on the size of the Skyrmions 
(related to d) as compared to the magnetic length I, and 
thus comparison with experiment and HF calculations 
will only make quantitative sense when c!> I. For "rea- 
sonable" parameters 
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these take the values I = 8.2 nm and d/l 7. 

The Skyrmions up to topological charge three have 
been studied in Ref. [ fl4"| , using radially symmetric 
Ansatze. In this paper it was shown that, for a large 
range of parameters, the charge-two Skyrmion is bound 
whereas the charge-three Skyrmion is not. They also 
made a comparison with more microscopic Hartree-Fock 
results, and discussed the limitations of the classical ap- 
proach. Several studies of the local baby Skyrme models 
show that the non-radially symmetric multi-Skyrmions 
can have far lower energies than their radially symmetric 
analogues. However, the exact details are highly depen- 
dent on the form of the potential term; see (T^j . Based on 
our past experience, we think it is very important to do a 
full two-dimensional numerical study. We first study the 
minimization of the energy functional (|^) . We discretize 




FIG. 1: The Skyrmions from topological charge one to three, 
for the parameter d/l = 7. Lighter colors denote a higher 
charge density, and lengths are expressed in units of d. 



the model on a rectangular spatial grid and impose fixed 
boundary conditions with the vacuum field n = (0, 0, 1) 
on the boundary. We start out with a field configura- 
tion of a given topological charge. Then we minimize 

We 



the energy functional using simulated annealing 1 1 



have extended this minimization technique to handle the 
long-range interaction terms arising from the Coulomb 
force. 

We have investigated the minimal energy solutions 
with the topological charge one to four; see Fig. |l| 
for some results. We found that only the charge-two 
Skyrmion is bound. The simulations for charge three 
and four seem to reveal a bound charge-three Skyrmion 
and a isosceles-triangular charge-four Skyrmion. How- 
ever, these are finite box effects. When increasing the 
box size we are left with individual Skyrmions of charge 
one and two. As energetically expected, a charge-three 
field splits into a charge-one and a charge-two Skyrmion, 
and the charge-four Skyrmion splits into two charge-two 
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FIG. 2: Magnetisation in triangular charge- two Skyrmion crystals, for increasing densities. The units of length is given in terms 
of the parameter d; for d/l = 7 the change in filling of the quantum-Hall state is (from left to right) Av — 0.024, 0.0544, and 
0.0967. Light colors denote a spin pointing upwards; the darkest denote a downwards spin. The contour levels are the same 
across all three graphs. See also |l6 . 



Skyrmions. We see clear indication that the charge- two 
Skyrmion might be the building block of any quantum 
Hall Skyrme crystal. It is very likely that a change in 
the parameters of the Lagrangian, and modification of 
the potential are able to dramatically change the struc- 
ture of the multi- Skyrmions and possibly lead to stable 
Skyrmions above charge two. In fact, any change of the 
effective action, such as higher-order gradient expansion 
terms, medium modifications to the Coulomb potential 
(see below), quantum effects, and inhomogeneities can 
have an impact on weakly bound multi-Skyrmions. 

We now turn our attention to the study of the quantum 
Hall Skyrme crystal, and produce a full two-dimensional 
study. We do not impose any symmetry, apart from 
that imposed by the periodic boundary conditions of 
the lattice. We need to be very careful when study- 
ing crystals with non-local Coulomb interaction terms. 
The crystal could be thought of as densely packed multi- 
Skyrmions, the quasi-particles of the Skyrme model, 
but the Skyrmions are all charged, and the samples 
used to study Skyrmions are electrically neutral. This 
means that Skyrme crystals exist in a positively charged 
background. This also renormalizes away the infinite 
Coulomb energy otherwise carried by a crystal. Since 
we shall use a finite size box for the simulation, we can 
use an Ewald sum rule to evaluate the sum over all the 
image charges generated by the periodic boundary con- 
ditions. This method explicitly gives us the zero-wave 
length component of the Coulomb field which we then 
renormalize away (a more detailed analysis of the proce- 
dure will be given in Ref. p2[). The first two parts of 
the energy expression (|5|) are not modified when apply- 
ing periodic BCs, but the Coulomb part becomes a sum 
over the vectors R — (n±L±,n2L2), where L\^ 2 are the 



lengths of the simulation volume, 
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Using Ewald's technique, we split the sum into a sum 
over the real and the reciprocal lattice (the integers k 
and I are used in the sum over the real lattice, k' and I' 
for the reciprocal one) 
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where 



rfc( v / s) 



The term k' = I' = is infinite. Adding a positive back- 
ground charge leaves us a constant energy shift Q . Us- 
ing (^) minus this shift is simple and straightforward. 
The sums all converge extremely rapidly. 

We have performed a study of these finite density quan- 
tum Hall systems, for the same parameters as the multi- 
Skyrmion simulations. Of course, we are already biased 
by our study of the multi-Skyrmions and expect a lattice 
of deformed charge-two Skyrmions. Our best guess is a 
triangular lattice, because it is the preferred lattice of the 
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local baby Skyrme model [|1TJ. Our goal was two-fold: to 
find the energetically favored symmetry of the lattice of 
charge-two Skyrmions and to show that other lattice con- 
struction are not stable. We started out with many dif- 
ferent configurations, mostly lattices made out of charge- 
one Skyrmions. Nearly all configuration minimized to a 
configuration made out of charge-two Skyrmions. Tak- 
ing finite box effects into account, the lattice seems to 
be triangular. Indeed we believe that the energetically 
preferred crystal is the triangular lattice made out of de- 
formed charge-two Skyrmions; see Fig. ^|. There appears 
to be a rather deep connection between these results and 
the complex (CP1) parametrization of the Skyrme field 
in terms of Weierstrass functions, which will not be pur- 
sued here, see however fl2| . 

The thermodynamics of the quantum Hall Skyrme 
model is another interesting challenge. What are the 
phases of the model at finite density and temperature? 
Our simulated annealing already contains the thermody- 
namic partition function, so it is well suited to tackle 
this question: see the phase portrait of the baby Skyrme 
model H]. However, the model differs in several impor- 
tant aspects. The baby Skyrme model has a local fourth- 
order term rather than the global Coulomb term. The 
fact that the dynamics of the Skyrme field is very dif- 
ferent may have an influence on the thermodynamics as 
well. There other kinetic terms have also been discussed, 
see Ref. Nonetheless, we expect some similarities, 

because the crystals are similar. 

The most thorough discussion of the phase diagram of 
Skyrmion crystals based on a "simple" potential Ansatz 
is given in Ref. pi] . Many interesting phases are dis- 
cussed, and we have tried to simulate a few of those, by 
starting a simulation in a periodic box with a Skyrme 
field of the right symmetry. In all cases we found a rapid 
decay to crystals of charge-two Skyrmions, or extremely 
frustrated crystals (due to the small box used). This does 
not mean that such crystals cannot occur, but it seems 
to be a generic feature of Skyrme-like theories to have 
a toroidal charge- two Skyrmion that is bound: It seems 
to be this feature that is driving our results. Thus we 
would like to argue that at sufficiently high densities the 
Skyrme crystal looks like the ones shown in Fig. ||. At 
lower densities the behavior is different. As long as we 
stick with the effective action (|l|) we shall see isolated 
clusters with charge two; if we modify the interaction, 
we expect connected networks, or even random clusters 
of multi-Skyrmions, as discussed in Ref. 111]]. We intend 
to investigate these questions further. 

The occurrence of a fluid at finite temperatures, which 
might well be the phase transition seen in experiment, 



is also of relevance. It has been argued that quantum 
fluctuations play a crucial role here as well fll9| , so the 
question about the nature of the fluid and the phase tran- 
sition deserves further study. 
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